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We study the Casimir interaction between a sphere and a cylinder both subjected to Dirichlet, 
Neumann or perfectly conducting boundary conditions. Generalizing the operator approach devel- 
oped by Wittman [IEEE Trans. Antennas Propag. 36, 1078 (1988)], we compute the scalar and 
vector translation matrices between a sphere and a cylinder, and thus write down explicitly the 
exact TGTG formula for the Casimir interaction energy. In the scalar case, the formula shows man- 
ifestly that the Casimir interaction force is attractive at all separations. Large separation leading 
• term of the Casimir interaction energy is computed directly from the exact formula. It is of order 

~ hcRi/[L'^ln{L/R2)], ~ hcRlRl/L^ and ~ ftc7??/[L* ln(L/i?2)] respectively for Dirichlet, Neu- 
, mann and perfectly conducting boundary conditions, where Ri and R2 are respectively the radii of 

■ the sphere and the cylinder, and L is the distance between their centers. 
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dl' I. INTRODUCTION 

> 

^ I Casimir effect [H has attracted a lot of attention due to its wide applications in different areas of pfiysics and its 
^ . potential impact to nanotechnology. For a good review about tfie subject, one can read for example tfie book 
Q^' Naively, tfie Casimir energy is a sum of the ground state energies of all the eigenmodes of the quantum field. However, 
the computation of the Casimir energy is not a simple task since the naive summation is divergent and regularization 
is required. Before the turn of the century, the exact computations were limited to some simple configurations with 
^ \ a particularly large amount of works being devoted to the configuration of two parallel plates. 

• Since the work of Lamoreaux [3] in 1997, Casimir force have been measured with high precision in various config- 
\ urations This has stimulated theorists to study the Casimir force between any two objects. At the beginning, 
, several approximation schemes were developed for this purpose, such as the semiclassical approach [1-0] , the optical 
' path approximation (sl-fioj and the multiple refiection approximation 11). 



Since the year 2006, exact computation of the Casimir energy between two objects has become possible. In |12h16| . 
Gies et al. derived a worldline representation of the Casimir interaction between two objects imposed with Dirichlet 
, boundary conditions. In ^17 ], Bulgac et al. computed the Casimir interaction energy between Dirichlet spheres or 
[ • between a Dirichlet sphere and a Dirichlet plate using the multiple scattering approach, whose application in Casimir 
J> ] effect can be dated back to the work of Balian and Duplantier [3 [3 • Later several other researchers computed the 
Casimir interaction in various configurations such as sphere- sphere, sphere-plane, cylinder-plane and cylinder-cylinder 
[20| - |3^ using some form of multiple scattering approach. A general scheme for computing the Casimir interaction 
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, between any two objects was developed by Emig et al. in 123] fo r scalar fields and in [26| for electromagnetic fields. 

A totally different approach was used by Dalvit et al. |35l - l37| . They used mode summation approach to compute 
the Casimir interaction energy of two eccentric cylinders with Dirichlet, Neumann or perfectly conducting boundary 
conditions. Recently, we developed a general scheme for computing the Casimir interaction energy between two 
objects from the perspective of mode summation approach [331, which is fundamentally equivalent to the scheme 
developed by Emig et al. from the perspective of multiple scattering theory. 

The basic ingredients in the exact representation of the Casimir interaction energy between two objects arc the 
T-matrices of each of the objects and the translation matrices between the objects. The T-matrix of an object 
only depends on the coordinate system chosen and the boundary conditions imposed on that object, and it has 
been calculated for a plane, a sphere and a cylinder under various boundary conditions. The more difficult part in 
obtaining the exact representation of the Casimir interaction energy is in the computation of the translation matrices. 
For cylinder-cylinder and sphere-sphere configurations, the translation matrices are well known. For cylinder-plane 
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and sphere-plane configurations, they have been obtained in [26|, |33|- To the best of our knowledge, the translation 
matrices for the sphere- cylinder configuration have not been worked out before, and hence the Casimir interaction 
between a sphere and a cylinder has never been studied. The goal of the current work is to fill in this gap. 

In Section |lTl we generalize the operator approach in [39'] to compute the translation matrices between a sphere 
and a cylinder for a scalar field, and hence write down the TGTG (exact) formula for the Casimir interaction energy 
between a sphere and a cylinder which are both imposed with Dirichlet boundary conditions or Neumann boundary 
conditions. In Section IIIIl we do the same for electromagnetic Casimir interaction between a perfectly conducting 
sphere and a perfectly conducting cylinder. In Section llVl we use the exact TGTG formulas to derive the leading term 
of the Casimir interaction energy when the separation between the sphere and the cylinder is lar ge. In the Appendix 
1X1 we use the proximity force approximation [i^, |4l| and the derivative expansion proposed by |43l | to compute the 
small separation leading order term and next to leading order term of the Casimir interaction. In principle, these 
can also be computed from the exact formula but it is too complicated to be included here. This problem will be 
addressed in the future. 



II. EXACT SCALAR CASIMIR INTERACTION ENERGY 




FIG. 1: The configuration of a sphere and a cylinder outside each other. 



In this section, we derive the exact scalar Casimir interaction energy between a sphere and a cylinder. As shown 
in Fig. [1] the centers of the sphere and the cylinder are situated at (0,0,0) and (L,0, 0) respectively. The radii of 
the sphere and the cylinder are given respectively by Ri and i?2. The length of the cylinder is denoted by H, and we 
assume that H ^ Ri, R2. 

In spherical coordinates {r^9^(j)) with center at O = (0,0,0), the scalar field (y5(x, t), x = {x,y,z), where x — 
r sin 9 cos 4>, y — r sin 9 sin (j), z — r cos 9, can be expanded as 
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Here 

c 



^z\^,k) = crh['\kr)YUd.'i^). cr = (1) 



where 

are the spherical Bessel functions, y;,„(6',0) are the spherical harmonics and P™{z) are the associated Legendre 
functions. The constants C\^^ and Cf"* are chosen so that 

On the other hand, the scalar field 'p(x, t) can also be expanded in cylindrical coordinates (p, </>, z) with center at 
O' = (L,0,0): 

^ E {a^k.Vnl{^,k) + h^u.vlt{^.k))e-^'^\ (2) 
where x = x' + O', x ~ x' + L — pcos0 + L, y = y' ^ psini/), z = z'. 

The constants C^°s and C°"* are chosen so that 

Using multiple scattering formalism [2^ or mode summation approach [38| . it was shown that the scalar Casimir 
interaction energy can be written in the form 



h 



OO 



Ecs / '^^Tr In (I - M{iO) , (3) 







where 

M(iO = Tl(^e)Ul'(^e)T2(^0U''(^e)• 

([3]) is called the TGTG formula. The matrices Ti and T2 are obtained by matching the boundary conditions on 
the sphere and on the cylinder respectively. These matrices are related to scattering matrices [U l23j | and have been 
computed in a lots of literature. They are diagonal matrices. For Dirichlet boundary conditions, the diagonal elements 
TL and r2^^ are given by 

1 (Rik) 



(Rik) 



In [R2^n^+kl 



Kn(^R2^/^^^Tkl 
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For Neumann boundary conditions, 

-i/,, 1 (Rik) + RikI', (Rik) 



-\Ki^i {Rik) + RinK'^^^ (Rin) ' 



Here 

c 

The hardest part of the problem in is to compute the translation matrices U^^ and U^^. This has not been worked 
out before. In the following, we generalize the operator approach introduced by Wittman [s^ (see also (SSJ) to 
calculate these translation matrices. 

As in [H, [3^ , we introduce the differential operator Vim defined by 



2l + l{l_rn)l / ^^y- („) 



47r {I + my. \ ik J ' \ik 

Here to > 0, Pi{z) is the Legendre polynomial of degree I and Pi"^\z) is its m-times derivative. It follows from the 
definition of spherical harmonics that applying to e*'' '", where k — k^Gx + ^y^y + k^ez, r — xe^ + yey + ze^, we 
have 

7',™e'''---r/„,(0fc,0fc)e^'^-'- (4) 
for / = 0, 1, 2, . . ., —l<m< I. On the other hand, one can show by induction that 

VhnMkr) = i'ji{kr)Yi,nie,<j)), Vimh^o\kr) = h[^\kr)Yim{e , ^) . (5) 

It can then be shown that the spherical wave functions ip\°^{:x., k) and ip°^{x, k) have the following integral represen- 
tations (see [ssl |39|): 

(^J^^f (x, k) =^ ^ d0fe ^ dOk smekYim{ekAk)e''''\ 

nout pco poc ik^x+ikyy±i^ -k'l -k^z (6) 

^Z\^,k)=:^ dkx dkyYimi9k,cPk) , , ^ ^ 0. 

For the cylindrical wave functions ip^^^ (x', k) and (x', fc), we introduce the operator Q„ as in [3l|. For n> 0, 

' dx' + idy' 



Qn = 
Q-n = 



dx' - idy' 

ikj_ 



It follows that 

for all n, where k.± = kxGx + kyBy, p = x'bx + y'ey. One can then show that (see [38|): 

nout poo ±i.y/ fc^ — feja;'+ifea y'+ik^, z (7) 

(^r^! (x, fc) =^ / dkye^-^- == , X ^ 0. 
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The translation matrices U"'^^ and U^^ are defined so that their elements and U^l satisfy the following 

identities: 



oo I 



1=0 m=-l 



(8) 



n— — oo 



To find Ul^^ , notice that (g]) and © imply that 



C 



rcg „2Tr 



(-1)" 



iPl"rn"'Pirn) (0) ^J—J J d(j)kj dOk Sm OkYhniOk , (t'k)Yl'' rn'' (Ok , (t>k) ^ S„,,, ^^„,C ^ 

Therefore, applying Vi.-m to both sides of the first equation of ^ and set x = 0, we have 



(9) 



12 ^ 1^ 



(-i)™(n-™^?.t)(-o',fc). 



Using ([7]) and (HJ, we then find that 



^k,0k e' 



1.2 L2 

A. 1 rLy 



Passing to imaginary frequency with k = in gives 

(/ - m)! ifcz 



(10) 



For U'^ let us consider the case where m > 0. The case where m < is analogous. Using ^ and ([7]), the 
second equation of (|5]) can be written as 



Cr'^-'Vimhi'^ (fc|x' + O' 

oo 



n— — oo 



27r '^«fe-''"27ri" 



0fc g*-^ k^ —k^ cos (f)f^x' k'-^ —k^ sin (p^y' -\-ikzz' 



In imaginary frequency, we have 



n— — oo 



(-1)" 



v!_ \ pi™-) 



/j^gin^fc g-y^K^+fef cos i^fci'- Y^K?+fcf sin (pky'+ik^z' 



K / kIx' + O' 



2 y 47r (/ + m)! V « 
The left hand side is an inverse Fourier transform. Taking Fourier transform, we find that 

r.27r 



' — ' -1 /'ZTT 



sin <t>^y 



7r(-l)" / 2^ + l(/-m)! r d^, + idy> 
' 2iJ Y 47r (/ + m)! 1^ ^ 



O' 
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Applying integration by parts, we find that 



Hence, 



dz' 



' ' K / kIx' + O'l 



dz'- 



-k|x'+0' 

kIx' + O'l 



_ p(™) 

=(-ir 



ik \ r°° p-Ky/(x'+L)^+y'^+z'^ 



K J K 



p,™ — i - 



ik \ I f°° g-\/'^^+'^'i+kl{x'+L)+ikyy' 



dk 



^J^^~+kf+k^ 



°° 1 p2-R 



ik, 

K 



dkii 



-yjK^+kl+kl{x'+L)+ikyy' 



.^n^ + kl + kl 

Applying the operator Q_„ to both sides and setting x' ^ y' — 0, we find that 



2iJK y 47r (/ + m)! ' V 



Unk^,lmi^C) — " 



«;2 + A;2 + ;;2 



7r(-l)" 2l + l{l-m)\ 



Hk 



ik. 



4. a + (^v/^). 



(11) 



Now we can write down the formula for the scalar Casimir interaction energy between a sphere and a cylinder. It 
is given by where the trace Tr is 



1=0 m=-l 



and 



n— — oo 



ikz 

K 



Making a change of variables 



we have 



k-, = Ksinhfl, 



/'OO 

^ j dO cosh 6IP;" (i sinh 61) is:„_™ (Lk cosh 9) T^^^ PJJ^' {-i sinh 9) K„r> (Lk cosh 0) 



(12) 
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The formula for Mim,Vm'{iC) contains P™(±isinh0) and P™'(±i sinh6') which are complex. This leads to 
some doubt whether the Casimir interaction energy ([3]) is real. In the following, we show that this is indeed the case. 
Moreover, the Casimir energy is always negative. 

Expanding the logarithm in ([3]), we have 



s=l ''^ i=l \h=Omi = -h/ i=l 



Ec.s - - 1^ dm\l^ III (13) 

with the convention that Is+i ~ h and m^+i = mi. 

Now for any and for any positive real number z, I,j{z), Il(z) and K^(z) are positive whereas Kl{z) is negative. 
Moreover, for ^ > 0, 

-^Ii+i{z) + zl^iz) > 0, + < 0. 

Hence, we find that for either Dirichlet boundary conditions or Neumann boundary conditions. 



Using the representation (|12l) . we find that 



n^'.™.,'.+i™.+iW) (14) 



can be written as 



n^'.™.,'.+i™.+i(*0 = n / positive function x J|P™'(^sinh0,)F;"'+^(-^sinh^?, 



Moreover, the positive function is also even in each of the variables 6*^. 
Using the formula (see [l^) for associated Legendre functions: 

priz) =i-irl^priz), 

we find that for m > 0, 



m > 0, 



Pr(±«sinh0) =(±*)'+'"-i-cosh™0 y (^) ,}^\ ^^'^\, sinh'-^^-"( 



2';! ^ Vj7 (' - 2? - m)! 

Notice that except for the factor (±i)'+™, the rest are real. Hence, in the product 

II P;7' {i sinh (?,)^r;r ^"^^ - 

the factor ih+mi-u+i-mi+i ^^y[ cancel off. This shows that (|T4)) is real. For integration over 9i, we will have a function 
that is positive and even in 9i multiply with 

(sinh6'0''-™'+''+i-'"'+^ (15) 

This function is positive and even if and only if {li — rrii + Z^+i — m^+i) is even, or in other words, {k — rrii) has the 
same parity as {U+i — m^+i). When {U — nii + k+i — m^+i) is odd, (|T5| is an odd function in 6i and therefore the 
integration over 9i vanishes. In summary, the function ()14p is positive if all (h ~ mi), i = 1, . . . ,s, have the same 
parity. Otherwise, it is zero. Hence, the representation ()13p shows that the Casimir interaction energy is always 
negative. 
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For the Casimir interaction force Fqhs, it is defined as the negative of the derivative of the Casimir interaction 
energy with respect to d, the distance between the sphere and the cyhnder. Namely, 

^Cas = — TTl-E'Cas- 

ad 

Notice that d only appears in L: 

L = Ri+R2 + d, 

and in (IT2|) . L only appears in the two terms Kn-m {Lk coshO) and K„i'^n {Lk cosh 9). In differentiating (|T4)) with 
respect to d, we will have 2s terms, each one is obtained by differentiating one of these modified Bessel functions 
which give a negative function. Combining with the same argument as before, we see that the Casimir interaction 
force is always negative. In other words, the force is attractive. This is true for all distances. 



III. EXACT ELECTROMAGNETIC CASIMIR INTERACTION ENERGY 



In this section, we consider electromagnetic Casimir interaction between a perfectly conducting sphere and a 
perfectly conducting cylinder. The electric field E and the magnetic field B can be expressed in terms of the vector 
potential A by 

E = -^, B = VxA. 

at 

In spherical coordinates centered at (0,0,0), the vector potential A{x,t) can be expanded as 

/oo ^ ^ 



Z— 1 m— — l 



where 



A™'*(x, k) =_4=V X ^L(x, k)rer 
V H' + 1) 

=C!fnkr)XUO,(t>), 

,TM,*/__ 

'k^l{l + l) 



A™'*(x, k) =-^===V X V X (p*„(x, k)r 



Here * = reg or out, with fl°^{z) — ji{z) and = h\^\z), 



are vector spherical harmonics (4^ . Obviously, 



Straightforward computation gives 



A™'*(x,fc) = ^VxA™'*(x,fc). (18) 



In cylindrical coordinates centered at (L,0,0), A(x,t) can be expanded as 

/oo poo 11 ^ 

E (a„..A™;-^(x',fc) + 6„fe,A™;™*(x',fc) 
-ca J -oo n=-oo 

+c„,,_ A™- -s(x', k) + d„,. A™- -*(x', k) 



(16) 



Ar™'*(x,fc)-TVxA™'*(x,fc). (17) 



(19) 
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where 



TE,* / / , \ -r-, * 



=C* ( —^9*n{k^p)ep - 5,* (fc±p)e0 ) e 



i7i(f)-\-ikzZ 



AL'!'*(x',fc)=^VxVx<,^e, 



A;, 



Here ^ and 5°"* (2) = ^^^(z). We have 



TM,*, 



,,/(x,fc)=-VxA™-(x,fc), 



(20) 



Using muhiple scattering approach [26] or mode summation approach |38| , one stiU finds that the Casimir interaction 
energy can be written in the form ([3]), but the components of the matrices T-'^, U^^, T^, U^-'^ are 2x2 matrices. For 
perfectly conducting sphere and cyhnder, the matrices and are well known. They are diagonal and each 2x2 
component is also diagonal, i.e., 



l.TE 



T, 





1,TM 
Im 



T, 



2.TE 

r 





2,TM 



with 



Ti2,TE/-t\ 



rp2,TM/ 



(Rik) + RikI[^^ {Rik) 
'^Ki^i (Rik) + RikK[^, (Rik) ' 



1.2 

rv ^ 



Kn[R2V^^+k^t 



The translation matrices U'^^ and U^"'^ have not been worked out before. The main results of this section are the 
explicit formulas for these translation matrices. We use the same approach as in the previous section, which is inspired 
by [3§|. First define the vector- valued operator [3§|: 



^ hn — 



1 



(21) 



where 



L = -r X V. 



In [39|, it has been shown that 

1 



yZ(ITT) 



2i 



(22) 



and 

A™' -«(x, k) = Cr^-^Vl^Mkr), A™' = Cr^-'VMkr) 

These imply that 



rticg p2TT pTT 

A™' '■°s(x, fc) =^ rfgfc sm0fcXi™(0fe, ^fc)e''' ■^ 

A™' ■■'=s(x, fc) =^ y^ d^k dOu sinOk'-^ X Xi„M, 



For the vector cyhndrical waves, consider the operator [s^: 

Qn = ^ (Qn+l - Qn-l) + ^ (Qn+l + Qn^l) 

2 At 

In [3I], it has been shown that 

A™;*(x',fc) = z-"C:Q„.gS(fc±P)^ 



It foUows that 



A™' -«(x', fc) =^ y^ (-e«J e-^'=e*-''+''=-, 

/^out roo ±iyj/k\^-k^x' +ikyy' +ik^z' 



/-"out /-oo ±i-i/ fe^ — fe^x'+ifcj,y' + ifc^z' 

A™-"*(x',A)=^ / *^(^e,Je»^^^ ^ , .^0. 

"'-00 Jk]_-kl 



The components of the translation matrices U^^ and U^^ are defined by 

-(X - O', fc) . 5: ^ «:r"A™^ ^^^(x, fc) + t/^:I^'™A™' -(X, k)) , 

00 I 

A™^ -(X - O', fc) . 5: «:r"A™' -(X, k) + Cr;™A™' -(X, k)) ; 



/ — I m— — l 



A™^ °-(x' + O', A.) = f i/ 1^ ^ (C/^1:?/"A™: -(x', fe) + C/^™™A™' -(x', fc) 



n— — 00 



Using (Hi]), (dH]) and ([201), it is easy to deduce that 

rrij.TE.TE _ rry,TM,TM ..ij,TE,TM _ ..ij.TM.TE ,-.--19nr91 
'-^lm,nk^ ~ '~'lm,nk^ ' ^im,n/c^ ~ ^im,nfe^ ' 'J — "-"^ ^J^- 

To find the translation matrix U^^, we use the fact that [ssj: 

QTCg „2ir j-TT f— 1)™+-'- 

3 

r2iT 



\ C f f ( — 1)™+ 

ri"m" ■ A™''"^j (0) =^ y^ dcj^k dOk sm9kXi„^^.iek,M ■ X/™(0fc, = ^^^ki"S; 
irr„,„ ■ A™'^"sj (0) =^ y^ dcj), dOk sin0fcXz--,„,,,(0fe, • f - x X^„(^^fe, <;ifc)j - 0. 



m" , — m 
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Apply the operator Vi.-m- to both sides of ([28)) and set x = 0, (l30l) imphes that 



12,TE,TE 



=|J(-ir+^(p,,-™-A™-') (-00 



27r(-l)'+ii 2l + l{l~m)\ 



+ V (/ + m) 



^fcy ^^Pr(cos0fc)ee. +isin0fcP/"'(cos0fc)e0. • (-ze^J e*("''")"^'= ^ 



1,2 _ 1,2 



rrl2,TE,TM 

lm,nkz 



1,2 _ 1,2 



Passing to imaginary frequency, we have 



..12,TE.TE..,x (-1)" /a^/o; , 1^([^lZ!^V^Z±Mp"' (t / .-2 i t.2 



V'I(ITT)V '(/ + m)! K 

.rl2,TE,TM,.^^ (-1)" /o, , -|^ (^""^)! "IK f ik^ \ ij / 2 ^ 1.2 



n— — oo 



(31) 



(32) 



(33) 



For the translation matrix U^^, we have to use the foUowings: 

, . ficg /•27r 

(Q„" • A™--^) (0) d<p, (-ze,J . (-*e,J e^(" +")^^ 

= _C-gr"<5„»,_„, (34) 
(q„" • A™^-^) (0) =^ d4>, i-^e,J ■ i-eej ^ 0. 

Notice that can be written as 

Cr'i-'Vi,nhi'^ (fc|x' + 0'|) 

E f°° dkz CJj^S f^'^ /Tr21.TE,TE / ■ \ , Tr21.TM,TE / n\ 

-n — , ^—00 ^0 



y gin^fc gi-y/ fc^ — cos k^—k^ sin (p^y' -{-ikzz' 

Let us concentrate on the case where m > 0. The case where m < can be worked out in the same way. As in the 
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scalar case, passing to imaginary frequency and taking Fourier transform, we find that 



— , -^0 



n— — oo 



X e e 



7r(-l)" / 2? + l {I -my. 
2Hk y 4ttI{1 + 1) {l + m)l\Y 



pm+l 



ik. 



Qrn+i + + - m + 1)P™-^ i Q„,_i 



ik. 



(35) 



2i 



pm+l 



iky 



— Qm+i -{1 + m)il -m+ 1)P™-^ Q„_i 



ik. 



ik. 



me.P,"' Q. 



+ kl + kl 



Here we have used the formula (22). Apply the operator Q_„- to both sides of ([55]) and set a;' = if = 0, ([M)) gives 

ik. 



tt21,TE,TE/ 



7r(-l)"+i / 2/ + 1 (;-m)! 
4i/K V 47r/(? + l) (/ + m) 



m+l 



+ (Z + m)(;-m + l)P;" 



K 



/ _|_ I ^.2 

7r(-l)"+i v/«^Tfcf / 2/ + 1 (I^m)! 
' K V 47r/(Z + 1) {l + my/ 



ik 



Applying 



^ (S-n+l + Q-n-l) + ^ (Q-„+l - Q-n-l) 

2 2i 



to both sides of ([35l) and set x' = y' = 0, one find that 



U. 



21,TM,TE 
nkz.lm 



7r(-l)"+ii / 2Z + 1 {l-m)\ n 
4Hk Y 47r?(/ + l) {l + m)\YA' ' 



P, 



m+l 



+ (z + m)(;-m + i)p; 



%k. 



_7r(-l) 



n+l 



+ fc2 y 47r?(/ + 1) (? + m)! 



K 



Hence, 



rT21,TE,TE/r ■/-\ 



{-ly+^TT / 2/ + 1 (/ - m)! + fc2 
mi \ 4ttI{1 + 1) {l + niy. ^ 



^in—n.kz 



p/M — ^ Wf„_„ LV^^Tfcf , 



(-l)"+i7r / 21 + 1 [l-my. mn 



-.PP 



K 



K,n-n [l^K^ + ki) 



(36) 



(37) 



(38) 



Hn y 47r/(? + l) (/ + m)! + fc2 ^ 
Finally, we can write down the formula for the electromagnetic Casimir interaction energy between a perfectly con- 
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ducting sphere and a perfectly conducting cylinder. It is given by ([3]) with 
T'"' j{2l + 1){21' + my. {V - m')\ 



Hm.V' 



2 V l(l + l)l'{V + 1) {l + m)\{l' + m')\ 



/ cosh ePl^'{i sinh e) r-^^r(* sinh 



^ nOQ I 

— —oo \ 



X 

n— — oo 



cosh^ 

777 

\ nnsh ft ^ ' 



X 



/ cosh 6iP,7"''(-i sinh 61) ---^P;7'' (-i sinh ( 
m 



cosh0 

m' 
cosh 9 




Kn-m {kL COShO) Km' -n {kL COSh 0) . 



\ —PP (~i sinh 6) cosh ftPP sinh ( 

^ cosh 6* 

The trace Tr in ([3]) is now given by 

oo I 
l—l m— — l 

where tr is the trace over 2x2 matrices. 

As in the previous section, one can show that the Casimir interaction energy is real. However, it is not straightfor- 
ward that it is always negative. 

IV. LARGE SEPARATION BEHAVIOR 

In this section, we compute the leading term of the Casimir energy when the separation between the sphere and 
the cylinder is much larger than their respective radii, namely, when 

L ^ Ri, i?2. 

We need to use the following leading behaviors for the modified Bessel functions: When z —?■ 0, 
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Making a change of variables 



\K,{z) + zK's{z) 3tt 
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1 fc poo / OO I2 \ s / . 



in p3|) . we have 

The large separation leading term comes from the term which has the lowest power of uj. For the case Dirichlet 
boundary conditions are imposed on both the sphere and the cylinder, the large separation leading term comes from 
the term with s = 1, /i = nii = 0. Namely, 

he [°° 

Moreover, the leading term of AI^q qo comes from the term with n = 0, i.e., 

Mnooo-l^^^T^ I de cosh eXo [uj cosh 6) f° ^^"^ ^"^^ Ko (uj cosh 6) . 
00,00 2Ki{au})J^^ ^ ^ is:o(6wcosh6') ^ ^ 



Here 



^1 , R2 



When i?2 ^ L, we obtain from ([55]) that 

/oo 
dfl cosh ^ii'o (w cosh 6') i^To cosh 9) . 
-00 

Therefore, the leading term of the Casimir interaction energy is 

HcR /"OO POO 

^Cas ~ - „ 9^9 1 /r / n ^ / ^ujuj / (19 cosh (w cosh 9) Ko {uj cosh ( 
27r^L^ln(L/i?2) Jo ^-oo 



47r2L2 1n(L/i?2) 7-00 cosh6' 

hcRi 
inL^ ln(L/i?2)' 

For the case Neumann boundary conditions are imposed on both the sphere and the cylinder, the large distance 
leading term comes from the term with s = 1, li — mi = or /i = 1, mi = 0, ±1. Namely, 

hr r°° 

Ec.s - " ^ / ^ Mo,00 + <0,10 + + ■ 

Now for each of these M^^ the leading term comes from n = 0, ±1. Using (|39| . we have 

poo r°° / \ 

/ / dcjwM d6'cosh3 6'(i^o(wcosh6')V2i^i(a;cosh6')^) 

Jo ' ^'^L^ Jo J-oo ^ ' 

_ \m\Rl 

45L5 ' 

poo Yf'Tf^ f°° f°° / \ 

/ duM^oio--^ dujLj^ d9 cosh^ 9 sinh^ 9 (Ko{ujcosh9f + 2Ki{ujcosh9f) 
Jo ' 47rL Jo J~oo ^ ' 

SR1R2 

" 15L5 ' 

/ dwMi^^i «^-4 / dujLj^ d9 cosh^ 9 (Koiujcosh9 f + Ki{ujcosh9 f + K2iu}COsh9 

Jo ' St^-^ Jo J-oo ^ 

llRlRl 
15L5 
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Hence, the leading term of the Casimir interaction energy is 



'^^'^ ^ 2^L \ 45L5 15L5 ^ ISL^ / 45^L6 ' 

For the case perfectly conducting boundary conditions are imposed on both the sphere and the cylinder, the large 
distance leading term comes from the term with s = 1, = 1 and mi = 0, ±1. Namely, 

he f°° 

i^cas '^^^'^ (M1040 + Mii,ii + Mi,_ia,_i) . 

For each of these M;„j^;„i, the leading term comes from n — Q. 

f°° 1 R^ f°° 

/ dwtrMioioW ^ — 7-^, — r/ dujuj^ d9 cosh^ 9 Kq (uj cosh 6 f 

Jo ' TT L-^ln{L/R2) Jo 

3L^ln{L/R2)' 

I rfwtrMi±ii±i« 5 — -J- — -/ dwwM d6'cosh6'(2sinh^6' - l)Xi (wcosh^)^ 

Jo ' 47rL'='ln(L/i?2) Jo J -00 

Rl 



12L^ ln(i/i?2)' 

Therefore, the leading term of the Casimir interaction energy is 

P he f i?! R^ 



" 2ttL y3L^ln{L/R2) ' 6LHn{L/R2) 
_ hcRl 
"~ 47ri4ln(i/i?2)' 

In summary, the large separation leading terms for Dirichlet, Neumann and perfectly boundary conditions are given 
respectively by 

poo.D hcRi 

^^"^ ^ 47rL2 1n(i/i?2)' 

N ^ _ Tlh^RfRl 
457rL6 ' 
hcRf 



jpo ^ 

^'"'^ 47rL4 1n(i/i?2)' 

It is interesting to compare these to the large separation leading terms between two spheres and between two cylinders. 
For two spheres with radii Ri and R2, the large separation leading terms for Dirichlet, Neumann and perfectly 
conducting boundary conditions are given respectively by [l3, [13, H^l : 

rnOO,D hcRiR2 



E. 



e: 



N mihcRlR^ 



leihcRfRl 

~ 967rL7 ' 
UihcRlR^ 



For two cylinders with radii Ri and R2 and length H, the large separation leading terms for Dirichlet, Neumann and 
perfectly conducting boundary conditions are given respectively by (2T| : 

poo,D fi^H 

- 87TL^\n{L/Ri)][\niL/R2)y 
poo,N _ _ 7hcHRlRl 



8^X2 [ln(L/i?i)] [ln(i/i?2)]' 
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For any boundary conditions, we find that the leading term of the Casimir interaction energy for the sphere-cyhnder 
configuration is intermediate between the sphere-sphere configuration and the cylinder-cylinder configuration. 

Numerical evaluations of the Casimir interaction energy Ecas for Dirichlet, Neumann and perfectly conducting 
boundary conditions are plotted respectively in Fig. [51 Fig. [3] and Fig. 21 and they are compared to the respective 
large separation leading term E^'^ derived above. The energies are normalized by Eq — hc/{2TTR) and they are 
plotted as functions of d/R for 1 < d/R < 100, where d is the distance between the sphere and the cylinder, and R is 
their common radius. In all the cases, we see that the large separation leading term agrees quite well with the exact 
value when d/R is large. 




Eo — hc/{2TTR) as a function of d/R. The dotted line shows the large separation leading term E^^^ normalized by Eq. The 
figure on the right shows the ratio of the exact Casimir interaction energy to the large separation leading term. 




V. CONCLUSION 



In this work, we discuss the Casimir interaction between a sphere and a cylinder, which has not been considered 
before. The sphere-cylinder configuration can play the same role as the sphere-plane configuration in Casimir exper- 
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boundary conditions normalized by Eo = hc/{2TiR) as a function of d/R. The dotted line shows the large separation leading 
term -E^l'g normalized by Eq. The figure on the right shows the ratio of the exact Casimir interaction energy to the large 

separation leading term. 

iments, since it does not have the problem of maintaining parallelism in the plane-plane configuration. We consider 
scalar interaction with Dirichlet or Neumann boundary conditions, and electromagnetic interaction with perfectly 
conducting boundary conditions. 

We first derive the exact formula for the Casimir interaction energy. The hardest part in the derivation is the 
computation of the translation matrices between the sphere and the cylinder. Motivated by the operator method used 
by Wittman [3§| to derive the translation matrices between two spheres, we have derived the translation matrices 
between a sphere and a plane and a cylinder and a plane in [ssl j. In this work, we use the same method to derive 
the translation matrices between a sphere and a cylinder. The results can be expressed in terms of modified Bessel 
functions and associated Legendre functions. The final exact formula for the Casimir interaction energy is quite 
complicated. Nevertheless, when both the sphere and the cylinder are imposed with Dirichlet boundary conditions or 
Neumann boundary conditions, we can deduce from the exact formula that the Casimir force is always attractive. 

Using the exact formula, we compute the leading terms of the Casimir interaction energy when the separa- 
tion between the sphere and the cylinder is large. It behaves like ^ fi,ci?i/[L^ ln(L/i?2)], ^ hcR\R\/L^ and 
~ hcR\ I [L"^ Ini^L / R2)\ respectively for Dirichlet, Neumann and perfectly conducting boundary conditions. Here 
Ri and R2 are the radii of the sphere and the cylinder respectively, and L is the distance between the centers of the 
sphere and the cylinder. As for other configurations, we find that at large separation, the Casimir force is strongest 
in the Dirichlet case, and is weakest in the Neumann case. 

Since Casimir force has been confirmed in the sphere-plane configuration, measuring Casimir force between a sphere 
and a cylinder should be possible. We believe that this configuration will play an important role in nanotechnology. 
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Appendix A: Small separation behavior 

In this section, we consider the small distance behavior of the Casimir interaction energy. In principle, one can 
compute the small distance asymptotic expansion of the Casimir interaction energy from the exact formula using the 
perturbative method developed in [IJ. However, this approach is very complicated and is beyond the scope of this 
work. Here we compute the leading order term from the proximity force approximation and the next to leading order 
term using the derivative expansion proposed in [4^ . 
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As is well-known, in the small distance limit, the leading term of the Casimir interaction energy is given by the 
proximity force approximation pollZH. Based on the cylinder, which we parametrize hy x = R2 cos (j> + L, y = R2 sin cj) 
and z — z, where < (j> < 27r and —H/2 < z < H/2, the proximity force approximation gives 



H 1—1 ^2 

dz d(t>£i^^m,z)), 

H J„_cos-i-52. 



where L = Ri + R2 + d is the distance between the sphere and the cylinder, (d) is the Casimir energy density 
between two parallel plates, which is given by 

for Dirichlet or Neumann boundary conditions, and 



h{(t),z) ^ y/{R2 COS <j) + Ly + (i?2 sin 0)2 + z^ - Ri = ^ Rj + + 2R2Lcos(j) + z'^ - Ri 

is the distance from the point {x,y,z) on the cylinder to the sphere. After some computations, we find that for 
Dirichlet or Neumann boundary conditions, the proximity force approximation gives 



ijPFA ^R,J (Al) 

^Cas lUOd^ '\ Ri+R2' 



For perfectly conducting boundary conditions, the leading term is twice of (jAip due to the two polarizations of 
photons. 

By letting the radius of the cylinder approaches infinity, i.e., i?2 — >■ 00, we obtain the sphere-plane configuration. 
Indeed, in the limit R2 ^ 00, (lAip gives 

TT^hc 



1440^2 



which is the small separation leading term of the Casimir interaction energy between a sphere and a plane. It is also 
interesting to compare (|Aip to the leading term of the Casimir interaction energy between two spheres of radii Ri 
and i?2, which is 

TT^Hc R1R2 



1440^2 i?i + i?2 



We see that for all three configurations (sphere-sphere, sphere-cylinder, sphere-plane), the Casimir interaction energy 
behaves like ~ 1/^2 when d <C 1. 

For the next to leading order term, it has been a subject of much interest in these recent few years. For the 
cylinder-plate configuration, the small separation next to leading order term was computed in [27j from the exact 
representation of the Casimir interaction energy using perturbation method with careful order counting. This method 
was then extended to the sphere-plate configuration [m-'J?], the cylinder-cylinder configuration [4^ and the sphere- 
sphere configuration [4^. There is no doubt that this method can be extended to the sphere-cylinder configuration 
considered here, although some tedious work is required. Less than two years ago, another approach was proposed 
to compute the small separation next-to-leading order term of the Casimir interaction energy. In [50| . Fosco et al 
performed derivative expansion on the path integral representation of the Casimir interaction energy and obtained 
an expression for the next to leading order term in terms of the height profile. They have extended their method in 
5lL [52! but so far their results can only be applied when one of the objects is planar. Inspired by [50|, Bimonte et al 
43j proposed that the Casimir interaction energy has a derivative expansion of the form 

EEI = £ cf2a;4as(^)(l + Pi{H)VH^ ■ VH^ + P2{H)VH2 ■ VH2 

+ py, {H)WHi ■ VH2 + P- {H)z ■ (ViJi X VH2) + . 

where E can be taken to be the z — plane parametrized by a; = ix,y), z ~ Hi{x) and z — H2{x) are the height 
profiles of the two objects with respect to S, and H = Hi — H2 is the height difference. The leading term of 
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is precisely the proximity force approximation. Using the invariance of the Casimir interaction energy with respect 
to tihing the reference plane S, it was found that 

P-{H) =0, 

PAH)=l (l-i/^l^^i§^)-/3.(i7)--/^.(F) 
=2 - MH) - P2{H). 

For Dirichlet, Nemxiann and perfectly conducting bomidary conditions, (3 = j3i = j32 is found to be a pure number 
that only depends on the boundary conditions, which is given by 
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The latter implies that up to the next-to- leading order term, the Casimir interaction energy of two perfectly conducting 
objects is equal to the sum of the Dirichlet and the Neumann Casimir interaction energies in the same geometry [43j . 
In the sphere-cylinder configuration that we consider here, we can let E to be the plane z = 0, let the sphere be 
+ + (z — Li)^ = R\ and the cylinder be + (z + L2Y — ^2- simplicity, we assume that Li > Ri > 0, 
L2 > R2 > 0, and let L = Li + L2- Notice that L = Ri + R2 + d, where d is the distance between the centers of the 
sphere and the cylinder. Then 



Hi =Li - JRi -x^~ y\ H2 = -L2 + JR'2 - x\ 



H =L- ^Rf - a;2 - y2 _ _ ^2^ 



i?2_^2_y2' ^ ^ 

Vi?i • ViJ2 = - 



9 

X 



^Rl^x^-y^^. 
After some computations, we find that 



, 1 _itRi I R2 d 3 2 , „2 

""h^ ~ d^ y R1 + R2 V RMRi + i?2) U ' ' 



1 VH2 . VH2 =31/ 0.0^ 



d\lR2{Ri+R2f 

Therefore, derivative expansion shows that up to the next to leading order term. 



pDE _ aTT^hcRi I R2 ( d 2 , p2 

^Cas - ^440^2 V i?i + i?2 V RMRi + R2) \8 ' ^ ' 

+Pd ( + , + —^^^ - (2 - 2p)d- 



Ri{Ri 4- -R2) R2{Ri + -R2) / -Ri + R2 



an^HcR, I R2 /i_5^_^^2/3-l)^+(/3_r)^ 



1440^2 \ R1+R2 \ 8R1+R2 Ri V 8 / i?2 
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where = = 1 and = 2. More precisely, we have 



^DE,D TT^hcRi / i?2 A 5 



cas 1440^2 \l Ri+R2\ 8Ri+R, 
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^DE.N TT^hcRi I R2 5 d 



E"t' = 7-^\ T^^^ 1 - ^ 



'^"'^ 1440^2 V i?i + i?2 V 8R1+R, 
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^DE p TT^hcRi / i?2 Z^-, 5 (i 
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24i?2/ 




40\ d 




^) R'i 




20\ d 
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24 " 


i?2 


7 ION 
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24 ~ 


i?2 



C''^'' 720^2 \l Ri+R2\ 8R1+R, 

It will be interesting to compare these results to the results computed directly from the exact formulas obtained in 
Section |lT] and Section IIIII This will be considered in a future work. 
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